The phase diagram of the hexagonal lattice quantum dimer model 

R. Moessner/ S. L. Sondhi^ and P. Chandra^ 
^Department of Physics, Jadwin Hall, Princeton University, Princeton, NJ 08544> USA 
^ NEC Research Institute, 4 Independence Way, Princeton NJ 08540, USA 

(February 1, 2008) 

We discuss the phase diagram of the quantum dimer model on the hexagonal (honeycomb) lattice. In 
addition to the columnar and staggered valence bond solids which have been discussed in previous 
work, we establish the existence of a plaquette valence bond solid. The transition between the 
plaquette and columnar phases at v/t = —0.2 ± 0.05 is argued to be first order. We note that 
this model should describe valence bond dominated phases of frustrated Heisenberg models on the 
hexagonal lattice and discuss its relation to recent exact diagonalisation work by J.B. Fouet et al. on 
the Ji — J2 model on the same lattice. Our results also shed light on the properties of the transverse 
field Ising antiferromagnet on the triangular lattice and the classical Ising antiferromagnet on the 
stacked triangular lattice, which are related to dimer models by duality. 



PACS numbers: 75.10.Jm, 05.50.+q 05.3Q.-d 
; I. INTRODUCTION 

Quantum dimer models (QDMs) were introduced to 
describe the physics of Heisenberg antiferromagnets in 
a regime dominated by valence bonds .u This regime is 
best realised in cases where the conventional Neel state 
is destabilised by quantum fluctuations or prohibited by 
frustration. The most prominent appearance of such 
dimer models has been in the context of the superconduc- 
tivity of the cuprates, where the QDM on thfij square lat- 
tice was introduced by Rokhsar and KiyelsonEJ to describe 
the physics of the slaort-range flavouifl of resonating va- 
lence bond physicsjj Subsequently, Read and Sachdev 

] showed that QDMs arise naturally in certain "extreme 
quantum limits" of generalizations of SU (2) magnets to 
SU{N) or Sp{N) symmetry with large Na 

More recently, an exact duality between QDMs in 
d = 2 and frustrated Ising models in a weak transverse 
field was explored by the present authors.|j This mapping 
connects the QDM not only to transverse field Ising mod- 

[ els but also to a class of ferromagnetically stacked frus- 
trated Ising magnets in dimension d = 2 + 1, which are 
of independent interestp^s they also have a range of ex- 
perimental realisations.El The connection between QDMs, 
short-ranged RVB physics and Ising gauge theories has 
also been discussed recently byrnFradkin and two of the 
present authors (RM and SLSj.Q 

[ From the work to date, it appears that QDM on 
the square lattice does not realise the disordered phase 

I envisaged in the short ranged RVB scenario of high- 
temperature superconductivity. Rather, it is ordered ev- 

' ery where except at a poipt, a possibility already noted 
in RK's original work,EI and fleshod|--OHt-|frjiin vari- 
ous viewpoints by different authors .siiy'ollillla'El While 
some evideace has been presented dissenting from this 
scenario as will become clear in the following, we 
believe there no longer is any real basis for doubting it. 

The physics of the square lattice QDM is closely con- 
nected to the physics of the hexagonal lattice QDM - 



both lattices exhibit critical classical dimer correlations 
which can be traced, via height representations, to their 
bipartite nature. (In contrast the non-bipartite triangu- 
lar lattice exhibits disorderpJ classical correlations and 
an RVB phase in its QDM.Eil) This connection was dis- 
cussed first by Read and Sachdev (see Ref. ^) and has 
more recently been discussed by us in the context of a 
study of frustrated transverse field Ising models.B 

Whereas the general structure of the phase diagram 
for both problems has been in place for some time, set- 
tling the detailed structure has turned out to be diffi- 
cult. Approaches based on mappings to height models, or 
Landau-Ginzburg- Wilson theory contain undetermined 
parameters upon which the detailed nature of the or- 
dering pattern depends. From a point of view of numer- 
ics, the problems of simulating the quantum problem has 
been restrictive in that the studies thus far have been Ujp-i 
ited to diagonalisations of systems of rather small sizeja'tll 
from which subtle difference in correlations have been dif- 
ficult to read off. 

In this paper, we solve this problem for the the QDM 
on the hexagonal (honeycomb) lattice and map out its 
phase diagram. This is done with the aid of the above- 
mentioned mapping to a classical stacked triangular Ising 
magnet. This mapping gives access to a number of ana- 
lytical results but, more importantly, allows efficient nu- 
merical simulations of systems much larger than the ones 
previously studied. 

Our central result is that the QDM on the hexagonal 
lattice has three phases, namely a staggered, a plaquette 
and a columnar valence bond solid (VBS). The transi- 
tions between these phases are a first order transition 
between the columnar and plaquette phases and a com- 
bination of a first order and continuous transition in the 
case of the plaquette and staggered phases; the latter in- 
volves fluctations in the ground state on one side of the 
transition but not on the other. 

Wc further discuss the implications of these results for 
the properties of stacked triangular Ising magnets with 
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nearest- neighbour interactions for which our algo- 

rithm allows us to avoid some numerical limitations en- 
countered in previous studies, and where the extension to 
the QDM provides considerable insight into the stability 
of the low temperature phase. This phase turns out to 
be one with three inequivalent sublattices, one of which 
is disordered, in accordance with the results of some, but 
not all, previous studies. 

We also review the connection of the QDM under study 
to frustrated Heisenberg antiferromagnets on the hexag- 
onal lattice. Such magnets are prime candidate for being 
described by the quantum dimer model, and it turns out 
that the Heisenberg model with competing interactions 
does indeed seem to reaJise the order present in two of 
the phases of the QDM.E3 

Turning to the QDM, its Hilbert space consists of hard- 
core dimer coverings of the hexagonal lattice. The Hamil- 
tonian acts on each hexagonal plaquette of the lattice. It 
contains two terms, a kinetic (T) and a potential {V) one. 
The former generates a plaquette resonance move by ro- 
tating a triplet of dimers pv 60° (see Fig. p, in analogy 
to the benzene resonance.cil The latter is diagonal in the 
dimer basis and simply counts the number of plaquettes 
able to resonate ('flippable plaquettes'). 

The Hamiltonian of the QDM can thus be repre- 
sented as a sum over plaquettes of the following plaquette 
Hamiltonian: 

Hqdm = -tf + vV 

= t + h.c) + V + \a){a\) . 

It has one free parameter, namely the ratio of the 
strength of the potential and kinetic terms, v/t. 

The structure of this paper is as follows. In Sect. 0, 
we discuss the phases which one might expect to en- 
counter in the model under consideration. Sect. ^ con- 
tains a summary of the methods used to establish the 
results that follow. The numerical results on the QDM 
are presented in Sect. from which the phase dia- 
gram (Sect. 0) follows. We then discuss implications for 
the study of magnets, namely triangular stacked Ising 
(Sect. 0) and 5 = 1/2 hexagonal Heisenber g (Se ct. [Vll ) 
models. We close with a conclusion in Sect. VIH. 



topological sector if one can be obtained from another by 
strictly local rearrangements of the dimers.cl) 




(a) (b) (c) 

FIG. 1. Dimer patterns on the hexagonal lattice: (a) stag- 
gered, (b) plaquette and (c) columnar. The marked links 
have a high probability of being occupied by a dimer in the 
respective phases. Note that in each case, there are only two 
inequivalent sets of links. A dimer plaquette move effected by 
T consists of rotating the three dimers surrounding a plaque- 
tte (like the one labelled with a plus) by 60°. 

As one decreases v through the ground state moves 
into another sector, which contains an exponentially large 
number of dimer configurations. The two candidate 
phases in this sector are depicted in Fig. |^b and c; these 
are the plaquette and columnar valence-bond solids, re- 
spectively. In fact, for v/t — > — cx), one can see that 
the ground state will be the columnar state, as this max- 
imises the number of flippable plaquettes favoured by the 
potential term. 

The point v/t = 1 is the RK point where each equal 
amplitude superposition over a winding number sec- 
tor is a ground-^itate. An analysis in terms of height 
representations^ shows that there is a diverging corre- 
lation length as one approaches this point from v < t 
and that the critical theory is Gaussian. In the same 
language the two candidate states mentioned above for 
v < t are flat but the competition between them cannot 
be settled in the same analysis. We now turn to an al- 
ternative mapping of the physics of the QDM which will 
allow us to settle that question by computation. 



III. USEFUL MAPPINGS AND NUMERICAL 
METHOD 



II. CANDIDATE PHASES 

As mentioned above, the QDM on the hexagonal lat- 
tice is closely connected to its square lattice version. 
Hence a number of known exact statements on the square 
lattice carry over mutatis mutandis to the hexagonal lat- 
tice. Firstly, for v > t, the ground state is the stag- 
gered state, \ip), depicted in Fig. |l|a. This follows from 
the fact that a lower bound on the energy per plaque- 
tte is mui{0,v — t}, and only \ip) saturates this bound 
for V > t, with Hqdm \^) = 0. The dimer configuration 
corresponding to \(p) turns out to constitute a topological 
sector of its own. (Two configurations belong to the same 



This alternative, duality, mapping crisply distinguishes 
between the different phases. This mapping takes the 
QDM in d = 2 onto a classical, stacked, frustrated, 
anisotropic Ising magnet in d = 2 + 1 on its dual lattice.El 
The Hamiltonian for that model reads: 

PcH Ising ^ K^'^ (JiCTj - if ^ a^ffi' + Pc'Vc XI ■ 
{ij) {ii') i 

Here, the a are Ising variable defined on the sites of 
a stacked triangular lattice; the sum on {ij) runs over 
nearest neighbour pairs in the plane, whereas the one on 
{ii') is over pairs in adjacent layers. Bi is the in-plane 
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exchange field experienced by spin i; if it is zero, the 
corresponding dimer plaquette is flippable. 

To generate equivalent Hilbert spaces, one has to take 
the limit of infinite exchange in the planes, — > +00, as 
there is a one-to-one correspondence between the hard- 
core dimer coverings on the hexagonal lattice and the 
Ising model ground states on the triangular lattice, up to 
a global spin reversal.o 

The equivalence then holds in the scaling limit K'^ 
-|-oo, with the quantum inverse temperature /3q given 
by fiqi = exp(2if'^)/7V = A/iV, where N is the num- 
ber of stacked layers, so that the zero temperature limit 
corresponds to a system with infinite extent in the stack- 
ing direction. The conversion of parameters between the 
classical (C) and quantum (Q) problem proceeds via the 
formula vq jt = (3cvc^- In the following, the quoted 
values oi v/t are to be understood as referring to the 
quantum problem. Note that A (which we will quote 
in the following) quantifies the discretisation error ~ it 
gives a rough measure of a typical correlation length in 
the stacking direction. 

For the case v = 0, this model has been studied in 
the past by several groups. In an. influential piece of 
work, Blankschtein and coworkeraij have carried out a 
Landau-Ginzburg- Wilson analysis for this model, which 
uncovered an Xl"-symmetric action with an (in d — 2 + 1 
dangerously irrelevant) six-fold clock term, which breaks 
the XY symmetry at sixth order. 

The ordering pattern obtained by this method is a 
three-sublattice ordering pattern; depending on the sign 
of the six-fold clock term, the three sublattices have the 
ordering pattern (-I-M, 0, — M), or (-|-m, —n, —n), where 
the amplitudes ^'J-wm and n are undetermined in the most 
general scenario.EZI Translating these back to dimer cor- 
relations, one finds that the former pattern corresponds 
to the plaquette VBS, whereas the latter corresponds to 
the columnar VBS. 

Simulating Hiding is evidently straightforward in prin- 
ciple using classical Monte Carlo simulations in d = 3. 
The only complication arises from the scaling limit which 
has to be taken, requiring long correlation lengths in tka 
stacking direction. Here, we use a cluster algorithm,E3 
the implementation of which is easy in Ising language, 
but would have been rather hard to guess at for simula- 
tions of a stacked dimer model. The mapping onto the 
stacked classical Ising model, together with the cluster 
algorithm, is what enables us to simulate system sizes 
which are substantially in excess of those treated so far 
in numerical studies of the QDM. 

IV. NUMERICAL RESULTS 

In this section, we display the results of our numerical 
simulations, which we argue demonstrate the existence 
of a phase transition from the columnar to the plaquette 
VBS around v/t = -0.2 ± 0.05. 



Since we know that the columnar phase is encountered 
in the limit v/t —f —00, we can rephrase the question 
of whether there also is the plaquette phase by asking 
whether there is a phase transition as v/t is increased 
towards -1-1. From the discussion in the previous section, 
it is apparent that, in spin language, the plaquette state 
has zero magnetisation whereas the columnar state does 
not. All we therefore have to do is to look for restoration 
of the Ising symmetry to discover whether there exists a 
plaquette phase. This we do in the following section. 

To give an impression of the general phase diagram, 
consider Fig. ^ where we plot the root-mean square mag- 
netisation, rrirms, of the equivalent Ising model, for small 
to moderate system sizes over a broad range of the pa- 
rameter v/t. We quote system sizes in terms of number 
of sites. This is the number of unit cells of the hexagonal 
lattice, which equals the number of dimers, and also the 
number of spins of the dual triangular lattice. From this 
plot, a number of important features are already visible. 
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FIG. 2. rrirms as a function of v/t, in the sector containing 
the columnar and plaquette phases. In this sector, the actual 
transition at ti/t = 1 to a state with strictly zero rrirms shows 
up as a depression thereof; Pgt = 0.083, A = 10. 

Firstly, for large —v/t, rrirms approaches its limiting 
value within the ground states of the Ising model, which 
is 1/3. Deviations are well-captured perturbatively in 
t/v, as depicted by the dashed line, which shows the low- 
est order result. Just to the left oi v/t = 0, one can 
witness the vanishing of rrirms, which gets sharper with 
increasing system size. 

Finally, at the Rokhsar-Kivelson point, v/t — 1, there 
is a kink in rrirms, which indicates the further transition 
to the staggered state. The staggered state has strictly 
zero magnetisation, but since it is in a different sector 
from the other two states, this shows up as a depression 
of rrirms as v is swept through t, as the system tries to 
accomodate local staggered correlations within the wrong 
topological sector as best it can. 



3 




FIG. 3. Enlarged view of phase transition region. Note the 
first order nature of the transition. The solid line is through 
the points for to the largest system size. fSqt = 0.083, A — 10. 

To get a handle on the details of the transition, we 
examine the same plot for a much narrower range of v/t 
near the vanishing of rrirms , for a wider range of system 
sizes, as depicted in Fig. |[ As the system size increases, 
the transition sharpens up into a discontinuous drop in 
rrirms- This drop separates the region on the left with 
rrirms decreasing with an almost constant slope, from 
that on the right, with rrirms being near constant and 
close to zero. The phase transition thus appears to be of 
first order, as will be discussed in more detail below. 

To underline this result, we plot the scaling of rrirms 
as a function of inverse linear system size for a number 
of values oi v/t near the transition in Fig. ^. This plot 
shows that rrirms settles down to a nonzero value on the 
left of the transition, at v/t — —0.25, whereas it scales 
to zero on the right, for v/t > —0.15. The transition is 
located around v/t = —0.2, where the scaling appears 
inconclusive. From this, we think it is conservative to 
estimate the transition point between the two phases to 
be located at v/t = -0.2 ± 0.05. 
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FIG. 4. Scaling of rrirms as a function of , the inverse 
of the linear system size. Pqt — 0.083, A — 10. 



We conclude this section by addressing potential sys- 
tematic errors arising from the introduction of the dis- 
cretisation in the stacking direction, since the mapping 
to the quantum dimer model is exact in the continuum 
limit only. 
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FIG. 5. Development of a function of A and 

I3q. The dashed line is for 5184 sites, the others are for 2304 
sites. Reducing the discretisation error (increasing A), and 
lowering the quantum temperature (increasing A'') sharpen 
up the transition. 

In Fig. ^ we show the plots of rrirms vs. v/t for a sys- 
tem of 2304 sites using different couplings in the stacking 
direction, K'^ , thus varying A, at a fixed quantum tem- 
perature. It can be seen that the transition sharpens up 
as A is increased, but moves only little as A changes from 
10 to 20. As the quantum temperature is lowered by a 
factor of two at A = 20, the transition sharpens further 
but again does not move significantly. These effects are 
therefore certainly within the error bars we give for the 
value of the critcal v/t. The case of the largest system we 
have studied (also displayed in Fig. ||) clearly also falls 
into this range. 

We note that the absence of finite-size effects at f = 0, 
upon increasing the number of layers, TV, at fixed (3c f^nd 
L, implies the existence of a gap in this part of the phase 
diagram. This is not surprising since at that point, we are 
far away from the phase transition, which is first order 
at any rate. However, this observation makes the exis- 
tence of a gapless excitation at this point, suggested in 
Ref. ^ seem rather unlikely. More generally, our results 
fit snugly into the expectations from the height represen- 
tation analysis as well the analysis of the transverse field 
Ising models (see below as well) and so there seems little 
doubt that the analysis in Ref. O is flawed. 
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V. THE PHASE DIAGRAM 

jhase diagram we have thus obtained is depicted 
The columnar-plaquette phase transition is of 
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first order, whereas the one at the RK point is a sec- 
ond order one, albeit with the somewhat peculiar feature 
that, coming from the right, it appears to be first order as 
no fluctuations are visible leading up to the critical point. 
However, coming from the left, a gap closes, giving rise 
to the gapless resonon excitations. I 



columnar 



plaquette I staggered 
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transverse field point RK point 
FIG. 6. Phase diagram of the quantum dimer model on 
the hexagonal lattice. The nature of the ordered phases is 
indicated above the axis. 

There are a number of theoretical reasons which lead us 
to conclude that the transition from plaquette to colum- 
nar VBS is first order, as the simulations suggest. Within, 
the framework of the Landau- Ginzburg- Wilson theory,E£l 
the critical point corresponds to the vanishing of the co- 
efficent of the six-fold clock term, so that the system 
could in principle fluctuate between all the degenerate 
XY states (including the columnar and plaquette ones) 
without encountering any barriers. However, higher 'har- 
monics' (clock terms) will presumably come into play as 
they are unlikely to vanish at exactly the same point as 
the leading one; it is these which will prevent the barriers 
between the plaquette and columnar state from vanish- 
ing. 

Further, we note that the symmetry groups of the two 
VBSs are not such that one of them is a subgroup of an- 
other, which would be a criterion within Landau theory 
for a continuous transition. This is in fact a somewhat 
subtle point as both phases break translational symme- 
try and retain a six-fold rotational symmetry. However, 
when trying to form domains of one phase within an- 
other, it turns out that the centres of rotational symme- 
try lie in distinct places for the two phases. 

This point, incidentally, is somewhat simpler in the 
square lattice, where the columnar phase breaks trans- 
lational symmetry in one direction and also rotational 
symmetry, whereas the plaquette phase breaks transla- 
tional symmetry in both directions but retains a four-fold 
rotational symmetry. 



VI. STACKED MAGNETS 

Our simulations apply equivalently to the hexagonal 
dimer model and to the stacked triangular magnets. We 
therefore briefly digress here to note some implications 
of our results to the latter system, on which a good deal 
of work has been done, in great part inspired by the ex- 
istence of experimental compounds realising this model; 
for a review of both theory and experiment, see Ref. ^. 

Recall that Hising at = reduces precisely to this 
model; there, the presence of the plaquette VBS cor- 



responds to the three-sublattice (-I-M, 0, — M) ordering 
pattern for the triangular magnet. This agrees with the 
results of Refs. p^-pO|, whereas it partially disagrees with 
Refs. |l5|,|| 

This result is somewhat surprising, as even in the low- 
temperature limit, it appears that one ends up with only 
a partially ordered state (spins on one sublattice are 
still equally likely to point up or down in this phase), 
rather than the apparently more fully ordered state 
(m, — n, — n). However, note that in either case, fluctu- 
ations are present down to zero temperature - in fact, 
these statea,|are stabilised by those fluctuations in the 
first place.tZlEl In the absence of fluctuations, the energy 
of any ferromagnetically stacked ground state configura- 
tion of the triangular magnet would be the same. Such 
a selection of an ordecjed state by fluctuations is known 
as order by disorder .E3 It has the feature that, although 
weak fluctuations arc needed for stabilising the state, 
their strengthening will lead to a melting of the order 
they themselves established in the first place. 

In the present case, it now so happens that the inter- 
mediate phase with a disordered sublattice can benefit 
from the fluctuations, and survive them. As fluctuations 
are suppressed (e.g. by adding a negative v, or a mag- 
netic field, see below), one enters the phase with a higher 
degree of ordering. 

We can be reasonably confident that upon lowering the 
temperature even further, there will not be a transition 
to the (m, — n, — n), mainly because the critical v/t seems 
to move very little, if at all in the right direction, as the 
temperature is lowered. Nonetheless, we are not entirely 
clear how to resolve the discrepanciesjyith Refs. |5 16. 
As for the hard-spin mean-field theoryjij it is conceivable 
that the fluctuations are somewhat underestimated there, 
thus landing it on the wrong side of the fine dividing line 
between the two states. At any rate, we have explicitly 
simulated temperatures lower than the expected transi- 
tion temperature, and found no transition: for a system 
of size 27 X 27 X 1024 spins at = 2.3 (A ~ 100), 
we find sublattice root-mean square magnetisations of 
(0.95|rXj.l2, 0.95). The early simulations by Blankschtein 
et aZJi3 may have run into problems as ergodicity is lost 
for single-spin updating at low temperature with an ip^ 
creasing correlation length in the stacking direction,L2l 
a problem only since remedied with the developmeal 
more advanced Monte Carlo simulation technology£ 
However, our results are fully consistent with the mean- 
field analysis of Ref. provided the coefRcicnt of the 
clock term remains of the same sign throughout the or- 
dered phase. rin 

To illustrate the closeness of the two phases,t3'E3 in 
Fig. we display the sublattice correlation matrix for 
a system of 2304 spins in a stack of height 120. It 
can be seen that upon application of a small magnetic 
field, one leaves the (-f Af, 0, —M) phase and enters the 
(to, — n, —n) one. Note that in the classical model, the 
proximity of the phases is artificially enhanced as the 
strength of the field is effectively multiplied by A as the 
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correlation length in the stacking direction increases. For 
this reason, the abscissa of the plot is scaled up by A. 

sublmag 48_1 20 




0.04 0.06 

FIG. 7. Correlation matrix of sublattice magnetisations 
(Mi,M2,M3), with Ml > M2 > Ms, as a function of mag- 
netic field, for a stack of 120 layers of 2304 sites. K"^ = 1.15, 
so that A = 10. Fox h = Q one finds the (+M, 0, -M) phase, 
which rapidly gives way to the (m, —n, — n) phase. 



VII. RELATION TO THE HEISENBERG 
ANTIFERROMAGNET 

The QDM was conceived to describe the physics of 
a S* = 1/2 Heisenberg antiferromagnet in a phase domi- 
nated by nearest neighbour singlet bonds (valence bonds, 
or dimers) . The question thus arises under what circum- 
stances a VB phase is energetically competitive compared 
to the Neel state. Note that a spin can form a valence 
bond with only one of its neighbours, whereas it gains en- 
ergy from all its antialigned neighbours in a Neel state. It 
is via the resonance moves, captured by T, that it might 
make up the energy difference. 

One step towards destabilising the the Neel state is to 
maximise quantum fluctuations, that is to say, consides 
spin S = \/2 systems. The original idea of AndcrsonO 
was to choose a frustrated (triangular) magnet, as this 
does not permit a Neel state in the first place. In addi- 
tion, in a lattice with low coordination the disadvantage 
of each spin forming only one VB is relatively less severe, 
thus favouring a VB state. 

Moreover, it is advantageous not to have closed loops 
which are very short. Tliis can be seen from within the 
framework of the QDMjd which, formally, is an expan- 
sion in the overlap between distinct VB configurations. 
This overlap is exponentially small in the number of VBs 
in which they differ. The shortest closed loop of even 
length on the lattice on which the VBs reside determines 
the lowest order term in this expansion. This length 
is four for the square and triangular lattice, but six for 
the hexagonal lattice, thus favouring the latter. Further- 
more, introducing frustrating further-nearest neighbour 



exchanges yields a smaller prefactor (but not a smaller 
expansion parameter) in the RK expansion. 

The hexagonal lattice with frustrating further neigh- 
bour exchangesE^ would thus seem to be a good candi- 
date for being described by some sort of QDM. In fact, 
this expectation appears to bephorne out by very re- 
cent work of J. B. Fouet et al.f^ who did exact diag- 
onalisation studies on a J1-J2-J3 S = 1/2 Heisenberg 
model on the hexagonal lattice. They indeed found va- 
lence bond phases. For frustrating J2/J1 = 0.4, they 
found a staggered phase, which gives way to an (at least 
short-range ordered) columnar or plaquette phase around 
J2/J1 = 0.3; the numerics, while not bepg entirely con- 
clusive, is suggestive of the latter phase .0 

This is in keeping with the QDM analysis, which in- 
deed suggests that this transition corresponds to cross- 
ing the RK point between the plaquette and staggered 
phases. We should note, however, that carrying out the 
perturbation theory within the dimer manifold in the 
spirit of Ref. |^ does not place the point v — t in be- 
tween those two values of J2/J1, and so a more micro- 
scopic prediction of the properties of Heisenberg models 
will probably require going beyond this by including any 
renormalizations needed to obtain an effective QDM with 
nearest neighbor bonds only. 



VIII. CONCLUSIONS 

We have mapped out the phase diagram of the QDM 
on the hexagonal lattice. We have established the exis- 
tence of a plaquette VBS intermediate between a colum- 
nar and a staggered one. This was achieved by combin- 
ing a duality mapping to a classical model with a Monte 
Carlo cluster algorithm. In the process, we were also 
able to go well beyond the limitations of previous numer- 
ical work on the stacked triangular magnet to confirm 
the nature of the low temperature phase in that model. 
Based on the close connections between the properties of 
dimer models on the hexagonal and square lattices ,eI we 
expect the square lattice QDM, which can be [Obtained 
from similar simulations on a stacked magnet,E3 to be- 
have in an analogous manner. As detailed in Sect. VII for 
the hexagonal lattice, this could shed some light on tke 
properties of a Ji — J2 square lattice Heisenberg model,L£l 
for which the derivation of the RK model is totally anal- 
ogous. Nevertheless the competition between the plaque- 
tte and columnar phases is a matter of microscopic detail 
and so a thorough study of the latter would appear to be 
in order definitively to settle the issue. We expect to 
tackle this question in future work. 
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